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The  experimental  measurement  and  modeling  of  liquid  chemical  agent  spread  and  sorption  on  a  porous  substrate  are 
described.  Experimental  results  with  the  nerve  agent  0-ethyl  5'-(2-diisopropylaminoethyl)  methylphosphonothiolate 
(VX)  demonstrate  that  the  wetted  imprint  volume  increases,  even  after  the  sessile  drop  volume  is  exhausted.  This 
indicates  the  wetted  imprint  is  only  partially  saturated,  and  a  multiphase  flow  problem  formulation  is  needed  to  predict 
the  VX  fate  in  porous  substrates.  Three  characteristics  and  their  changes  in  time:  (i)  sessile  volume  remaining,  (ii)  wetted 
imprint  area  on  the  sand  surface  where  the  droplet  is  deposited,  and  (hi)  VX  penetration  depth  into  sand,  are  computed 
numerically  and  compared  to  experimentally  measured  values.  A  very  good  qualitative  and  quantitative  agreement  was 
found  between  the  numerical  and  experimental  results.  These  numerical  and  experimental  methods  can  be  used  to 
determine  the  spread  and  sorption  of  hazardous  materials  into  a  variety  of  substrates. 


1.  Introduction 

The  interactions  between  a  liquid  and  a  substrate  are  an 
important  factor  in  determining  the  overall  fate  of  the  liquid.  In 
the  case  of  a  liquid  deposited  onto  porous  substrates,  the  droplets 
sorb  into  the  substrate,  from  where  they  may  evaporate,  react,  or 
be  sequestered.^”^  Previous  studies  of  the  sessile  droplet  fate  have 
mainly  addressed  the  sorption  time  and  spread  during  periods 
when  there  is  sessile  volume  present  at  the  porous  medium 
surface.  The  sorption  time  and  the  shapes  of  the  sessile-drop  vol¬ 
ume,  and  the  wetted  imprint  have  been  predicted  using  analytical 
models"^’^  and  numerical  simulations.^’^  In  these  models,  the 
capillary  pressure  at  the  free  interface  is  obtained  from  the 
Washburn  equation,  and  the  wetted  imprint  saturation  is  set 
equal  to  one  (i.e.,  fully  saturated).^  The  liquid/substrate  interac¬ 
tions  influence  the  sessile  volume  shape  and  contact  angle,^’^^ 
where  the  sorption  dynamics  may  be  determined  by  the  liquid  slip 
velocity.^ The  wetting  influences  the  capillary  force,  which  is 
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opposed  by  the  drag  force  of  the  porous  substrate,  with  the 
resultant  of  those  forces  governing  the  sorption  process. 

However,  once  there  is  no  sessile-drop  volume  remaining,  the 
hquid  can  continue  to  spread  into  previously  dry  regions  within  the 
porous  substrate,  thus  increasing  the  wetted  volume.  Therefore,  the 
wetted  volume  of  porous  substrate  is  not  fully  saturated  by  liquid, 
i.e.,  the  pores  are  not  filled  with  liquid.  Instead,  the  liquid  proceeds 
along  preferential  flow  paths  producing  a  partially  wetted  imprint 
of  volume  much  larger  than  the  fully  saturated  volume.  Along  these 
hnes,  the  sorption  times  for  some  pairs  of  liquid/powders  can  not  be 
predicted  using  the  fully  saturated  imprint  assumption,  indicating 
the  formation  of  a  partially  saturated  imprint.  The  experimental 
evidence  from  MRI  measurements  further  corroborate  this  par¬ 
tially  saturated  imprint  multiphase  flow.^^”^^  Using  the  capillary 
network  model,  Markicevic  and  Navaz^^  have  shown  numerically 
the  transition  between  fully  and  partially  saturated  imprint 
during  the  droplet  sorption.  This  implies  that  the  droplet  sorption 
should  be  solved  as  a  multiphase  flow  problem,  with  the  multi¬ 
phase  flow  parameters,  capillary  pressure  and  relative  perme- 
ability^®”^^  included  in  the  momentum  conservation  equation.^^ 
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Figure  1.  Schematic  of  the  sessile  droplet  sorption,  including  the 
definitions  of  inlet  Ai^x  and  interface  Ax^t  boundary,  dynamic  6^, 
and  liquidi solid  6  contact  angle  for  sessile  volume  Vses-  (a)  greater 
than  and  (b)  equal  to  zero. 

The  determination  of  multiphase  flow  parameters:  capillary 
pressure  and  relative  permeability  have  been  described  by  Navaz 
et  al.  for  the  blister  agent  bis(2-chloroethyl)  sulfide  and  agent 
simulants  on  a  range  of  porous  substrates.^^ 

The  experiments  in  this  article  were  carried  out  with  the  nerve 
warfare  agent,  0-ethyl  *S-(2-diisopropylaminoethyl)  methylphos- 
phonothiolate  (VX),  sorbing  into  medium-grain  sand.  The  sessile- 
drop  volume,  wetted  imprint  surface  area,  and  penetration 
depth  were  measured  in  time.  The  time  in  which  the  experimental 
results  were  collected  was  much  longer  than  the  time  it  takes  for 
sessile  drop  to  sorb  into  the  sand.  In  this  way,  the  numerical 
and  experimental  results  are  compared  over  a  large  number  of 
experimental  points  rather  than  for  one  experimental  point  (i.e., 
the  time  it  takes  for  liquid  to  infiltrate  a  porous  medium).  The 
multiphase  sorption  influences  other  transport  rates  such  as 
evaporation,  reaction,  and  the  agent  fate  significantly.^^’^^  Given 
the  large  number  of  liquid/substrate  combinations,  it  is  useful  to 
develop  a  predictive  computational  model  to  obtain  insight  into 
this  process.  Such  a  model  would  be  applicable  to  pesticides, 
industrial  chemicals  and  have  environmental,  defense,  and  home¬ 
land  security  implications. 


2.  Model  System 


angle  within  the  porous  medium,  6  defines  the  angle  of  wetted  and 
nonwetted-regions  and  remains  the  same  throughout  the  sorption 
process. 

The  liquid  flow  in  the  porous  substrate  is  formulated  as  a 
multiphase  flow  problem  for  gas  and  liquid  phases.  The  problem 
is  described  by  coupling  continuity  equation,  Darcy’s  Law,  the 
capillary  pressure  jump  condition  at  the  gas/liquid  interface,^^ 
and  the  sum  of  phase  saturations  (^cas)  and  (^uq): 


=  0  (1) 

u''>  =  -pl'lg-x)  (2) 

Pc  =  PGas  -  PLiq  (3) 

^Gas  ^Liq  ~  1  (4) 


where  [/]  can  be  either  the  gas  or  liquid  phase  of  viscosity  /t  and 
phase  permeability  K^.  The  phase  velocity  u  and  pressure  p  are 
obtained  from  eqs  1-4  with  the  capillary  pressure  Pc  equal  to  the 
pressure  difference  of  nonwetting  and  wetting  phases.  The  in¬ 
stantaneous  position  of  the  free  interface  area  Ax^^x  is  determined 
from  the  liquid  phase  velocity.  The  change  of  the  wetted  volume 
Fwet  in  time  is  calculated  after  integrating  the  liquid  phase  flow 
over  Ax^{. 


wet 


[Liq] 


(5) 


An  initial  sessile  drop  volume  Lses(^o)  of  liquid  is  deposited  onto 
the  porous  medium  surface,  where  in  the  sorption  process,  the 
instantaneous  volume  remaining  on  the  surface  decreases, 
and  sorbed  volume  Lsor(0  increases  in  time.  The  sorbed  liquid 
wets  the  porous  substrate  of  porosity  0,  where  0  is  the  ratio  of 
void-space  volume  to  total  volume.  The  wetted  imprint  volume 
is  equal  to  V^et  =  ^sor/0  for  fnlly  saturated  sorption,  or  L^et  > 
Lsor/0  for  nonsaturated  sorption.  The  spread  and  sorption  are 
driven  by  a  capillary  force,  and  even  when  there  is  no  liquid  left  on 
the  porous  medium  surface  VgesiO  =  0,  the  sorption  continues  and 
(Fwet)  increases  in  time  with  the  volume  saturated  by  the  liquid 
decreasing  at  the  same  time.  The  capillary  force  is  responsible  for 
the  shapes  of  the  sessile-drop  and  the  wetted  imprint,  and  a 
coupled  problem  of  flow  in  the  surface  liquid  together  with  the 
flow  in  porous  substrate  has  to  be  solved.  In  some  special  cases, 
such  as  constant  droplet  base  radius,  the  flow  problem  solution  in 
the  sessile  volume  is  known,  and  only  the  flow  within  the  porous 
substrate  needs  to  be  determined.  Figure  1  shows  the  liquid  fate 
during  the  spread  and  sorption  into  a  porous  substrate  for  the 
special  case  of  the  constant  droplet  base  radius  (rq).  The  region 
where  liquid  sorbs  into  the  medium  is  defined  as  the  inlet 
boundary  area,  Axni.  The  region  where  the  liquid  front  advances 
through  the  substrate  is  defined  as  the  liquid  interface  area,  Axnt. 
The  dynamic  contact  angle  changes  in  time  and  is  defined  as  an 
angle  between  porous  medium  surface  and  tangent  at  the  sessile 
droplet  shape  from  the  contact  line.  This  differs  from  the  contact 
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Except  for  the  surface  of  substrate  on  which  the  droplet  is 
deposited,  all  other  boundaries  are  defined  as  no-flow  for  either 
phase.  This  condition  transforms  into  a  zero-pressure  gradient  in 
the  direction  normal  to  the  boundary.  For  the  boundary  on  which 
the  sessile  volume  resides,  the  gas  phase  flows  out  from  the  porous 
substrate  domain  and  through  the  boundary  fraction  that  is  not 
covered  by  liquid.  It  is  important  to  note  that  only  a  part  of  this 
boundary  is  covered  by  liquid,  defined  by  A^^x  in  Figure  1 .  The  gas 
pressure  is  equal  to  the  atmospheric  pressure  p^xm^  and  there  is  no 
liquid  flow  across  this  surface.  At  A^^,  the  liquid  pressure  is  equal 
to  the  pressure  within  the  sessile  droplet,  and  there  is  no  gas  flow 
across  the  area  defined  by  A^^.  Once  there  is  no  sessile  liquid  left 
^ses  =  O5  gas  pressure  becomes  equal  to  />atm  over  all  domain 
boundaries,  and  liquid  pressure  at  the  free  interface  is  calculated 
from  eq  3. 

The  capillary  pressure  in  eq  3  is  calculated  from  the  Leverett 
function,  scaled  by  substrate  porosity  0,  single-phase  flow  perme¬ 
ability  Xsp,  and  contact  angle  0.  The  phase  permeability  of 
each  phase  is  expressed  from  and  the  relative  permeability 
that  is  a  quadratic  function  of  saturation.  The  details  of  the 
parameter  determination  and  numerical  method  based  on  the 
continuum  approach  is  described  in  Navaz  et  al.^"^  In  brief,  the 
governing  equations  describing  mass  and  momentum  conserva¬ 
tion  (eqs  1  and  2)  for  liquid  and  gas  phases  are  written  over  a  small 
volume  and  then  assembled  and  solved  for  overall  computational 
domain.  The  partial  differential  equations  in  a  continuum  domain 
are  solved  on  a  finite-difference  mesh  with  fourth-order  Runge— 
Kutta  explicit  numerical  scheme  used  for  time  integration.  The 
algorithm  is  second-order  in  space  and  fourth-order  in  time. 
The  numerical  solution  includes  each  phase  saturation,  pressure 
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Figure  2.  Change  of  the  dynamic  contact  angle  6^  of  VX  sessile 
droplet  on  a  medium-grain  sand.  Inset  image  shows  the  sessile 
droplet  for  time  t  =  73  ms  after  freefall  impact.  A  red  outline  is 
provided  for  image  clarity,  and  the  width  of  the  white  box  repre¬ 
sents  2  mm. 

and  velocity,  and  how  each  changes  in  time.  From  the  solution, 
other  time-dependent  characteristics  can  be  obtained,  including 
wetted  volume,  area  of  imprint,  and  protrusion  depth  into 
substrate.  Initially,  the  numerical  solution  is  validated  against 
experimental  results  by  comparing  the  time  it  takes  for  volume  to 
infiltrate  porous  substrate,  and  a  very  good  agreement  was 
documented. 


3.  Experimental  Procedures 

Caution!  The  handling  of  chemical  warfare  agents  should  only 
be  performed  by  trained  personnel  at  an  approved  facility  using 
applicable  safety,  security,  and  surety  precautions. 

A  medium  grain-size  sand  AFS-50,  98%  Si02  acquired  from 
Warmwell  quarry  (Bardon  Aggregates)  with  the  following  prop¬ 
erties:  density  of  2.65  g/cm^,  BET  surface  area  of  0.33  m^/g,  and 
grain  size  between  0.25  mm  and  0.50  mm  was  used  in  experiments. 
The  samples  of  nerve  warfare  agent  O-ethyl  S-(2-diisopropyl- 
aminoethyl)  methylphosphonothiolate  (VX)  were  taken  from 
existing  stocks  and  were  determined  to  have  a  purity  of  at  least 
90.3%.  To  each  stock  aliquot  of  VX,  3  wt  %/wt  iisopropyl- 
carbodiimide  (Fluka)  stabilizer  was  added  to  simulate  a  muni- 
tions-grade  VX  stock.  In  all  experiments,  a  droplet  of  initial 
volume  equal  to  6  mm^  was  deposited  onto  sand.  In  the  set  of 
experiments,  the  spread  and  sorption  time  characteristics,  includ¬ 
ing  (i)  dynamic  contact  angle  0^,  (ii)  imprint  surface  wetted  area 
and  (iii)  penetration  depth  were  measured,  as  described 
below. 

For  the  dynamic  contact  angle  measurements,  the  VX  was 
deposited  as  a  6  mm^  droplet  using  a  micropipet  (Hamilton  10 //L 
adjustable  pipet,  model  1708—32)  allowing  the  droplet  to  freefall, 
rather  than  to  be  ejected  forcefully.  Since  the  infiltration  rate  of 
the  sessile  volume  into  sand  is  fast,  a  high-speed  video  camera 
(MotionXtra  HG-TH,  Redlake,  1000  frames/s)  was  used  to 
capture  freefall  velocity,  contact  angle,  and  sessile  droplet  lifetime 
on  the  sand.  Contact  angle  measurements,  0^,  and  droplet  volume 
calculations  were  obtained  using  custom  Matlab  routines  from  a 
series  of  side-on  sessile  volume  images,  as  shown  in  Figure  2. 

Wetted  area  and  penetration  depth  measurements  were  carried 
out  using  the  apparatus  shown  in  Figure  3,  both  schematic  (left) 
and  actual  design  (right).  The  sand  container  top  surface  has 
interior  dimension  of  9.2  cm  x  3.3  cm,  enabling  the  deposition  of 
multiple  drops.  The  surface  spread  area  A^^  was  captured  using  a 
video  system  (Sony  XC-ST50,  Imperx  capture  card,  custom  soft¬ 
ware  by  Battenkill  Technologies)  at  a  rate  of  once  per  3  min.  The 
penetration  depth  images  were  captured  using  a  digital  SLR 
camera  (Canon  digital  Rebel  EOS  SLR,  Canon  control  software). 


with  frames  captured  once  per  minute.  Images  were  analyzed  for 
the  area,  eccentricity,  and  equivalent  diameter  using  a  custom- 
written  Matlab  routine.  The  code  used  a  combination  of  auto¬ 
thresholding  and  morphological  processing  routines  to  identify 
the  droplet.  The  construction  of  the  sand  container  enables  the 
placement  of  the  camera  underneath  the  apparatus  such  that  it 
does  not  interfere  with  the  sorption  process. 

Once  identified,  the  geometric  parameters  were  calculated. 
Movie  frames  were  rendered  of  the  original  results  with  an  overlay 
of  the  perimeter  boundary,  to  enable  visual  comparison  between 
the  original  image  and  calculated  boundary.  An  example  of  a 
processed  movie  frame  for  A^  is  shown  in  Eigure  3  (top-down 
image).  To  facilitate  determining  d^,  the  construction  of  the  sand 
container  with  inclined  bottom  wall  creates  a  substrate  of  con¬ 
tinuously  variable  and  known  depths.  The  top-down  camera 
measures  the  spread  dynamics  of  the  drop  on  the  surface  A^, 
while  the  bottom- up  camera  determines  penetration  time  by  color 
change  of  the  M8  chemical  warfare  agent  detector  paper,  placed  at 
the  bottom  of  the  incline  before  adding  the  sand.  The  depth  of 
penetration  d^  is  calculated  from  the  known  slope  of  the  bottom 
wall,  a  =  7°,  and  the  distance  of  the  initial  color  change  measured 
from  the  box  leading  edge.  Examples  of  top-down  and  bottom-up 
images  (right)  used  to  determine  A^  and  d^  are  included  in 
Eigure  3.  Einally,  for  the  droplets  for  which  VX  did  not  reach 
the  indicator  paper,  the  wetted  imprints  were  taken  out  from  the 
sand,  and  from  their  images,  the  imprint  shape  was  calculated. 

4.  Results  and  Discussion 

Erom  the  experimental  results,  the  time  characteristics  of  the 
dynamic  contact  angle  6^,  imprint  wetted  area  A^,  and  penetra¬ 
tion  depth  d^  are  calculated.  The  change  of  in  time  is  calculated 
for  one  droplet  only.  The  characteristic  A^{t)  and  a  single  value  of 
d^  were  measured  for  each  droplet  applied.  The  experimental 
results  are  compared  with  the  numerical  predictions.  In  the 
numerical  solution,  the  density,  dynamic  viscosity,  and  surface 
tension  ofVX  are  set  p=  1.01  g/cm^,//  =  1.23  x  10~^  Pa  x  s,  and 
a  =  3.2  X  10~^  Pa  x  m.  The  medium  grain-size  sand  porosity 
and  permeability  are  determined  to  be  equal  to  0  =  0.44  and  K= 
6.05  X  10“^^  m^.^^  The  values  of  the  capillary  pressure  and 
relative  permeability  are  given  by  Navaz  et  al.^"^ 

The  droplet  fell  a  distance  of  2.3  mm  within  18  ms,  and  was 
observed  to  remain  as  spherical  during  the  descent  using  the  high¬ 
speed  imagery.  Because  of  impact  on  the  surface,  oscillations  were 
observed  within  the  droplet  for  the  initial  30  ms  after  impact. 
Hence,  the  dynamic  contact  angle  is  measured  starting  at  40  ms 
after  the  impact.  The  Weber  number  is  defined  as  the  relative  ratio 
of  kinetic  to  surface  energy.  This  number  can  be  calculated  from 
the  impact  velocity  (around  Wim  ^0.2  m/s)  and  the  characteristic 
length  (the  diameter  of  free  falling  droplet,  4ir  =  2.254  mm) 
according  to  We  =  pUu^d^^jo  which  is  around  We  3.5.  This 
relatively  low  number  indicates  that  the  droplet  impact  velocity 
did  not  overcome  the  surface  energy,  a  finding  further  supported 
by  the  nondistorted  spherical  shape  of  the  free-falling  droplet.  The 
sessile-drop  volume  after  40  ms  was  calculated  to  be  5.9  mm^, 
implying  that  no  significant  volume  sorbed  into  the  sand  in  this 
initial  time  interval  compared  to  the  initial  volume  of  6.0  mm^. 
The  volume  is  calculated  from  the  height  and  footprint  diameter 
of  the  sessile  drop  with  V=  l/3jr/z^(^/2<i  —  h),  with  a  typical  error 
of  ±10%.  In  Eigure  2,  the  results  of  measured  6^  are  shown  for 
time  up  to  t  =  1 92  ms.  Measurements  were  not  taken  once  =  25° 
as  sand  grains  are  comparable  in  size  with  the  sessile  droplet 
height  that  would  produce  large  errors  in  ^d  calculated.  Eor 
shorter  times  and  sessile  droplet  heights  being  larger,  the  sessile 
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Figure  3.  Experimental  setup  and  imprint  shape  measurements:  principle  of  penetration  depth  apparatus  (left)  and  an  actual  apparatus 
photograph  (right).  The  images  of  wetted  imprint  area  (top  right)  and  breakthrough  point  were  captured  by  indicator  paper  (bottom  right). 
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Figure  4.  Summary  of  wetted  imprint  surface  area  results:  (a)  typical  frame  and  its  overlaid  contour  (red  line)  used  to  measure  A^, 
(b)  individual  experimental  results  A^{t),  (c)  numerical  solution,  and  (d)  experimental  averaged  A^{t)  and  confidence  interval  (±2crAw)  given 
by  full  symbols  and  error  bars,  and  numerically  calculated  A^(t)  for  1%  (empty  squares)  and  2%  (empty  triangles)  saturation  threshold. 


droplet  shape  can  be  clearly  distinguished  as  can  be  seen  from  the 
inset  in  Figure  2,  for  t  =  73  ms.  In  Figure  2,  is  determined  using 
two  methods:  (i)  direct  measurement  of  the  angle  between  the 
substrate  and  the  droplet  outer  shape  finding  the  tangent  line,  and 
(ii)  calculated  from  the  ratio  of  the  drop  height  h  to  radius  of  the 
base  To,  where  tan(0d/2)  =  hjr^?^  The  agreement  between  two 
methods  is  excellent  (circles  and  triangles  in  Figure  2),  where  from 
both  methods  it  is  found  that  there  is  a  linear  decrease  in  the 
contact  angle  of  0.4°/ms.  The  ratio  method  is  appropriate  only  for 
sessile  droplet  in  the  form  of  a  spherical  cap.  The  observed 
agreement  between  two  methods  implies  that  VX  on  medium 
grain-size  sand  indeed  forms  spherical  cap  geometry.  The  sessile 
droplet  sorption  time  is  also  predicted  numerically,  where  the 
influence  of  deposition  velocity  is  neglected  and  the  droplet 
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base  radius  tq  remains  constant.  The  numerically  calculated  value, 
625  ms,  compares  excellent  with  the  experimentally  measured 
sorption  time,  613  ms. 

Once  the  sessile  VX  has  sorbed  into  the  sand,  a  further  wetting 
of  sand  by  VX  occurs  as  the  agent  sorbes  into  previously 
non  wetted  volume.  Thus,  both  A,^  and  increase  further  in 
time.  The  experimental  and  numerical  trials  are  carried  out  up  to 
250  min.  The  results  reveal  that  the  imprint  wetted  area  ^4^  ranges 
from  20  mm^  to  1 10  mm^,  and  penetration  depth  from  1  mm  to 
4  mm.  Thus,  an  increase  more  than  four  times  for  both  A^  and 
dp  is  observed,  giving  VX  saturation  ^  reduction  for  over  an  order 
of  magnitude  (s  <  0.1).  The  droplet  experiment  was  repeated  ten 
times  obtaining  ten  distinct  ^w(0  curves,  from  where  the  results 
repeatability  is  analyzed.  From  these  experiments,  the  average 
wetted  area  after  t  =  200  min  is  89  mm^,  with  a  standard  devia¬ 
tion  of  13  mm^  and  a  relative  standard  deviation  of  15%. 
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Figure  5.  Power  law  comparison  for  two  numerical  (^th  =  1  %  and 
5th  =  2%)  and  the  experimental  wetted  area  results.  The  changes  in 
powers  n  show  that  the  threshold  saturation  is  in  the  range  of 
between  1%  and  2%. 

The  variation  in  was  also  observed.  It  should  be  noted  here  that 
there  are  two  sources  of  the  variability  of  the  experimental  data. 
One  is  from  the  experimental  error  and  second  one  is  from  the 
sand  structure.  Because  of  the  droplet  size  and  pore  size  distribu¬ 
tion  in  the  sand,  depending  on  where  the  droplet  is  placed,  the 
liquid  surface  flow  (large  pores  beneath  the  droplet),  or  liquid 
protrusion  flow  (small  pores)  can  be  favored.  This  is  corroborated 
by  experimental  results  for  where  after  t  =  1  min,  in  two 
different  experiments  d^  was  found  equal  to  0.7  mm  and  1.1  mm, 
respectively. 

Figure  4  summarizes  both  experimental  and  numerical  results 
for  A^.  The  wetted  imprint  area  is  determined  from  experimental 
images  (part  a),  used  to  calculate  the  wetted  area.  The  individual 
experimental  results  are  shown  in  Figure  4b.  In  all  cases,  A^^ 
increases  rapidly,  and  then  levels  off  as  A^  approached  the 
maximum  spread.  The  individual  experiments  were  self-consistent 
and  smooth  compared  to  the  test-to-test  variance.  The  numerical 
results  are  shown  in  Figures  4c  and  4d.  The  colors  in  Figure  4c 
indicate  the  level  of  saturation  for  VX  in  the  sand  as  a  function  of 
spatial  position.  The  saturation  is  defined  as  a  fraction  of  pore 
volume  occupied  by  liquid  VX.  The  average  A^^  is  shown  in 
Figure  4d  with  full  red  circles,  and  bars  represent  the  95% 
confidence  interval  (ib2(7Aw)-  It  should  be  mentioned  here  that 
the  experimentally  measured  obtained  from  image  processing, 
does  not  provide  information  on  the  saturation  level  of  the 
medium,  nor  is  the  threshold  value  known  for  differentiating 
wet  versus  dry.  However,  in  the  numerical  solution,  this  satura¬ 
tion  threshold  for  the  wetted  surface  visibility  was  set  at  1  %  and 
2%,  with  A^{i)  profiles  represented  by  square  and  triangle 
symbols  in  Figure  4,  respectively.  The  numerical  results  indicate 
that  setting  the  saturation  threshold  smaller  gives  a  larger  wetted 
area,  with  both  numerical  solutions  (1%  and  2%  threshold 
saturation)  within  confidence  interval  of  A^(t)  computed  from 
the  experimental  results. 

The  numerical  results  from  Figure  4  show  that  the  change  of 
the  wetted  area  depends  on  the  threshold  saturation  (^th)  which  is 
defined  as  a  minimum  saturation  for  which  the  sand  becomes 
visibly  wet  (recordable  by  camera).  The  experimental  results  are 
between  the  1%  and  2%  saturation  threshold  curves  that  are 
obtained  by  numerical  solution.  For  all  three  profiles  A^(t),  the 
wetted  area  initially  increases  more  rapidly  compared  to  later 
times.  However,  for  the  longer  spread  time,  the  experimental  and 
one  of  the  numerical  results  (^th  =  2%)  for  the  wetted  area  reach 
almost  a  constant  value.  It  takes  a  longer  time  for  the  1% 
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Figure  6.  Changes  of  penetration  depth  as  a  function  of  time  dp(t) 
including  experimental  values  (full  circles),  numerical  predictions 
(empty  triangles),  and  square  root  dp(t)  dependence  (dashed  line). 
The  shape  of  wetted  imprint  and  length  scale  (5  mm)  are  shown  in 
the  inset. 

saturation  threshold  wetted  area  curve  to  demonstrate  the  same 
behavior.  In  order  to  investigate  these  changes,  the  wetted  area 
results  are  analyzed  using  a  power  law: 
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A^(t)- =  at"  (6) 

where  ^w,o  is  the  wetted  area  right  after  the  sessile  volume 
becomes  equal  to  zero,  and  a  and  n  are  the  power  law  constants. 
The  results  for  A^(t)  as  obtained  from  the  numerical  solutions 
(5th  =  1  %  and  5th  =  2%)  and  experiments  are  shown  in  a  log-log 
plot  in  Figure  6,  using  different  symbols.  The  symbols  are 
connected  with  thin  lines  for  better  visibility.  The  thick  solid  lines 
represent  the  power  law  linear  regression  from  eq  6.  The  powers 

^  0.47  and  ^  0.39  are  determined  for  the  experimental  and  the 
numerical  results  of  2%  saturation  threshold  wetted  area,  respec¬ 
tively.  Setting  the  saturation  threshold  to  5th  =  1  %,  the  numerical 
A^(t)  (square  symbols)  follows  the  power  law  in  the  whole  time 
range  with  «  ^  0.52.  This  implies  that  the  threshold  saturation  is 
chosen  properly,  and  its  value  is  in  the  range  1%  <  5th  <  2%. 
From  the  same  numerical  solution  as  shown  in  Figure  4b,  the 
penetration  depth  dp(t)  is  also  calculated. 

The  comparison  between  experimental  and  numerical  results  is 
given  in  Figure  5  with  circular  and  triangular  symbols,  respec¬ 
tively.  The  numerical  values  of  dp(t)  are  computed  for  saturation 
threshold  equal  to  1%.  It  can  be  seen  that  the  numerical  results 
closely  follow  the  experimental  data.  The  numerical  solution 
curve  dp(t)  is  not  smooth,  caused  by  a  computational  grid.  Since 
there  is  smaller  number  of  computational  nodes  in  the  penetration 
depth  direction,  the  nodes  with  saturation  less  than  the  saturation 
threshold  are  truncated  and  some  small  “jumps”  in  the  dp(t) 
profile  are  observed.  In  calculating  the  A^(t)  profile,  there  are  a 
larger  number  of  nodes  and  smooth  numerical  A^(t)  is  calculated. 
For  the  experimental  cases  in  which  VX  does  not  reach  the 
indicator  paper  placed  at  the  sand  apparatus  bottom,  the  wetted 
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imprint  is  removed  from  the  sand  and  photographed  as  an 
illustration  of  the  wetted-imprint  shape.  The  image  of  one  of 
those  imprints  is  shown  in  the  inset  in  Figure  5.  Additional 
investigation  of  the  imprint  shape  reveals  that  the  imprint  is  a 
spheroid  in  shape;  the  result  that  is  observed  from  the  numerical 
solution  also.  However,  it  is  important  to  note  that  the  is  not 
measured  from  these  wetted  imprint  shapes,  as  there  may  be  liquid 
present  outside  the  region,  but  not  enough  for  induce  clumping. 
The  protrusion  depth  is  modeled  as  a  power  function  of  time, 
where  it  is  found  that  is  a  time  square  root  function,  represented 
by  dashed  line  in  Figure  5.  Finally,  the  square  root  dependence 
(ip  of  time  t  has  been  reported^ for  the  dynamics  of  the  capillary 
driven  flows  into  porous  substrates. 

5.  Conclusions 

The  experimental  method  to  measure  the  dynamics  of  the 
chemical  warfare  agent  VX  sorption  into  medium  grain-size  sand 
for  both  sessile  droplet  present  and  after  no  liquid  is  left  on  the 
porous  medium  surface  is  developed.  The  sessile  volume  sorption 
rate  and  shape  are  determined  using  the  high-speed  video  camera, 
where  the  sessile  droplet  is  in  the  shape  of  spherical  cap,  and  the 
dynamic  contact  angle  decreases  in  time  linearly.  The  wetted 
imprint  is  a  spheroid  in  shape  with  its  base  referred  as  wetted 
imprint  area  and  other  axis  defined  as  penetration  depth.  Both  of 
these  quantities  are  measured  experimentally  by  recording  the 
surface  wetted  area  and  identifying  the  arrival  time  of  VX  to  the 
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sand  bottom  of  known  depth.  The  imprint  volume  increases  faster 
initially  compared  to  the  longer  times,  and  the  penetration  depth 
can  be  expressed  as  a  square  root  function  of  time.  The  infiltration 
process  is  also  solved  numerically,  in  which  besides  liquid  and 
porous  medium  physical  properties,  the  multiphase  flow  trans¬ 
port  parameters  capillary  pressure  and  phase  permeability  need 
to  be  known.  All  three  characteristics:  (i)  dynamic  contact  angle, 
(ii)  imprint  wetted  area,  and  (hi)  penetration  depth  are  computed 
and  compared  with  experimental  results  with  a  very  good  agree¬ 
ment  found  over  the  whole  time  interval  during  the  sorption 
occurs.  These  numerical  and  experimental  methods  can  be  used  to 
determine  the  spread  and  sorption  of  hazardous  materials  into  a 
variety  of  substrates.  Such  studies  are  applicable  to  pesticides, 
industrial  chemicals  and  have  environmental,  defense,  and  home¬ 
land  security  implications. 
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